Reheating via a generalized non-minimal coupling of curvature to matter 
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In this work one shows that a generalized non-minimal coupling between geometry and matter is 
compatible with Starobinsky inflation and leads to a successful process of preheating, a reheating 
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term for a scalar field \, and also encompasses a previously studied preheating channel based upon 
a non-standard kinetic term. 
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I. INTRODUCTION 

Cosmology is in a new era in which it is possible to 
make detailed quantitative analysis for the early uni- 
verse, due to the wealth of observational data stemming 
from important experiments such as WMAP [l[ and the 
Sloan Digital Sky Survey 0, amongst others. Inflation, 
which assumes a period of accelerated expansion after 
the Big Bang, is the most studied candidate to solve the 
monopole, horizon and planarity problems 0,0] , as well 
as explaining the anisotropics in the cosmic microwave 
background radiation. The simplest models of inflation 
yield a scale-invariant spectra of gravitational waves and 
energy density perturbations which act as seeds for struc- 
ture formation. Most models of inflation are based on 
the slow-roll of scalar fields. In contrast with the cosmo- 
logical constant scenario, slow-roll does not present the 
so-called "exit problem" , since the period of cosmic accel- 
eration is followed by a subsequent radiation-dominated 
era, along with a transient matter-dominated phase. 

Given the known equivalence between scalar models 
and the so-called f(R) theories, where the Einstein- 
Hilbert action is generalized to admit a non- linear f(R) 
term in the scalar curvature R, it is natural that inflation 
may also be obtained within the framework of the latter 
proposals (see Ref. Q for a thorough review) . In particu- 
lar, an early model for inflation relies on a quadratic addi- 
tion to the linear curvature term, f{R) = R + R 2 / (6M 2 ) 
[3] , with WMAP normalization of the CMB temperature 
anisotropics indicating that M ~ 3 x 10~ 6 Mp, Mp being 
the Planck mass ml. 
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A recent generalization of f(R) extensions of General 
Relativity (GR) includes the presence of non-trivial terms 
in R by introducing a non-minimal coupling of the scalar 
curvature with matter Q, via the action 



S = 



fi(R) + MR)C 



(1) 



This model has a rich lore of implications, both theoret- 
ical and observational. These include the deviation from 
geodesic motion Q, the possibility of mimicking dark 
matter by leading to the flattening of the galaxy rotation 
curves Q and the modelling of the accelerated expansion 
of the Universe at late times Although not related 
to this work, recent developments in inflationary models 
have led to a rehabilitation of the Higgs boson as a pu- 
tative candidate for the inflaton field, provided that it is 
non-minimally coupled to the curvature fioj . 

In the context of inflation, a non-minimal coupling is 
the key ingredient for the so-called preheating mecha- 
nism, which arises due to a scalar field x with a variable 
mass term of the form m 2 ^ = m + due to the oscil- 
latory phase that follows the slow-roll regime, the dynam- 
ics of this scalar field can undergo parametric resonance, 
thus giving rise to the quantum production of massive 
particles even for relatively low values of the coupling 
£, ^ 1 [HI El- Furthermore, preheating may also be im- 
plemented with more evolved couplings, e.g. a quadratic 
coupling R 2 X (HELl or via a non-standard kinetic term 
for the matter scalar field, of the form g(R)(dx) 2 [l5| . 

Hence, it is quite natural to expect that the model in 
Eq. ([I]) is capable of generalizing the preheating scenario 
to the case of a universal non-minimal coupling between 
matter and geometry. Thus, the main purpose of this 
work is to show that preheating indeed occurs if one as- 
sumes a generalized coupling of the form /2 {R)C in the 
modified Einstcin-Hilbcrt action, instead of the ad-hoc 
terms £i?x 2 or g(R)(dx) 2 - 
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However, the universality of this model implies that 
curvature is coupled to all matter species, i.e. matter 
and radiation, besides the scalar field x- Naturally, this 
may potentially modify the cosmological dynamics, as 
given by the Friedmann equation: the latter will display 
terms arising not only from the Starobinsky prescription 
fi(R) = R + R 2 / (6M 2 ), but also from the non-minimal 
coupling function f2(R)- Hence, one should first establish 
the validity of the Starobinsky inflationary regime, that 
is, that these extra terms do not become dominant during 
the inflationary phase. 

This manuscript is organized as follows: firstly, one 
briefly discusses the main features of the Starobinsky in- 
flationary regime and ensuing preheating mechanism, as 
well as the fundamental results arising from the consid- 
ered non-minimally coupled model. The third section 
addresses the required constraints to the non-minimal 
coupling that allow for Starobinsky inflation to occur. 
One then proceeds to the main purpose of this work and 
study how the non-minimal coupling enables a paramet- 
ric resonance leading to preheating of the Universe. Con- 
clusions are then put forward, and an outlook of future 
developments is drawn. 

Since it is not directly related to the obtained results, 
the analogy between the non-minimally coupled model 
under scrutiny and a multi-scalar-tensor theory [l(| is 
deferred to an appendix. A second appendix deals with 
the possibility of implementing an inflationary era purely 
via a non-minimal coupling between curvature and mat- 
ter. 



II. STAROBINSKY INFLATION 

As already stated, the well known Starobinsky infla- 
tionary model d| considers a quadratic correction to the 
Einstein-Hilbcrt Lagrangian, 



f(R) = 2k[R + 



R 2 
6M 2 



(2) 



where k = (WnG)- 1 = M^/Wtt and M ~ 3 x 1CT 6 M P . 
This model can lead to an inflationary stage in the early 
Universe, due to the presence of the quadratic term, 
which ends by the dominance of the linear one. 

One assumes the Robertson- Walker metric, as given 
by the line element 



ds 2 



-dt 2 + a 2 (t) 



dr 2 



\J\ — kr' 



dQ< 



(3) 



(4) 



Tpv + A M „Fi + — F\R)g llv 



with F\ = f[(R) and = V M V„ — g^JH- Introducing 
the Starobinsky prescription Eq. leads to 



R 
3AP- 



R 



^g^R ) = 



(5) 



1 

2k' 



i 



;A„ V R- 



R 2 



3M 2 12Af 2 
Inserting the metric Eq. (j3|), these may be rewritten as 



H 



IP 
2H 



1 



M 2 H 



-3HH 



R + 3HR + AI 2 R = 



(6) 
(7) 



where H = d/a is the Hubble parameter, R = 6(H+2H 2 ) 
and the dots denote time derivatives. 

The cosmic acceleration can be written as d/a = 
H 2 (l — e), where e = —H/H 2 is the slow-roll parame- 
ter. In the slow-roll approximation e <C 1, the former is 
positive defined, corresponding to an inflationary regime. 
Since in this approximation one also has \H / (HH)\ <C 1, 
the two first terms in Eq. © can be neglected, thus 
yielding 



M 2 
6H 2 



(8) 



One may straightforwardly calculate the number of c- 
folds N, so that the scale factor increases by an amount 
e N during inflation. This is given by 



= I Hdt 



2ei (ti 



(9) 



where t = ti and t — tf correspond to the instants of 
onset and end of inflation, respectively. 

Inflation ends when the slow-roll parameter Eq. © 
becomes of order unity, as the Hubble parameter drops 
below Hf ~ M/a/6. After this, one can no longer dis- 
regard the second time derivative of the scalar curvature 
in Eq. ([7]). For the solution of this equation one must 
perform the substitution R — > a~ 3 / 2 R, 



where a(t) is the scale factor and k is the spatial curvature 
(which is set to k = in the following calculations). 

The Einstein field equations of standard f{R) theo- 
ries are obtained from action Eq. ([T]), by substituting 
fi(R) = f(R) and taking the trivial case f2(R) = 1: 



R 



M 2 - -H 2 - -H ) 

4 2 



(10) 



During reheating one has M 2 » {H 2 , \H\}, which leads 
to a harmonic oscillator solution for the last equation 
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and a damped one for the original scalar curvature, 
R ex a~ 3 / 2 sin(Mf). In order to seek for the Hubble 
parameter one may neglect the r.h.s. of Eq. ([6]) ob- 
taining the solution H(t) oc cos 2 (M</2), which suggests 
the Ansatz 



H{t) = f(t) cos 2 (Aft/2) 



(11) 



Taking the solution for the slow-roll regime, H = —AI 2 /6 
and the approximation R(t) ~ 6H, one has (l7j 

~ 6£T = -3Mf(t) Bin[M (i - i„)] , (12) 
where i Q denotes the onset of the oscillatory regime and 



/(*) 



(13) 



In the M(t — t a ) ^> 1 regime one can consider the 
approximation 



/(*) 



3(t - 1„) ' 
for this last set of equations. This yields 



(14) 



H(t) 



3(t - 1 ) 



M 



(t - t„) 



(15) 



so that (H) ~ (2/3)(£ — 1 ) 1 (indicating that during the 
Starobinsky reheating phase, the Universe behaves as if 
matter-dominated) and 



R(t) ~ - 



AM 

t-t 



sm[M(t -t )} 



(16) 



A. Standard preheating 

Almost all the matter that constitutes the Universe, 
in the subsequent radiation-dominated era, was created 
during the reheating process at the end of the inflation 
through the gravitational particle creation, which occurs 
via oscillations of the Ricci scalar. 

This particle production can be effectively described 
by scalar fields 1 . For simplicity, one considers a scalar 
field x with mass m, along with a non-minimal coupling 
with the scalar curvature, 



1 For a discussion involving two scalar fields see e.g. Ref. [lS 



s = 



—gd x x 

f(R) i 



(17) 



2k 2 



2 2 

-m x 



where f(R) is equivalent to the function fi(R) in the 
action of our model. The variation of the action with 
respect to the field x leads to 



□X - m x -£Rx = 



(18) 



This quantum field x can t> e decomposed in modes as 
depicted below 



X(*,x) 



1 



(2tt) 3 / 2 



d 6 k x 



(19) 



a kX k(t)e- lk -* + a{xUty- 



where one has the creation a\. and annihilation ajt oper- 
ators of particles with mass m and momentum k. 

The field x can b e quantized in curved spacetime by 
generalizing the basic formalism of quantum field theory 
in flat spacetime. Then each Fourier mode Xk(t) obeys 
the following equation of motion 



(20) 



Xk + ZH X k + -» + m z + £R )Xk = , 



where k = |k| is the comoving wavenumber. After the 
substitution Xk a~ 3 ^ 2 Xk, this last equation becomes 



Xk 



^+ m 2 + £R- ^H 2 - ) 



The last two terms in Eq. ([21]) can be neglected if 
|£| > 1. Taking the approximation of the scalar curvature 
given by Eq. (|TB)) in the regime M(t — t a ) » 1, the same 
equation becomes 



Xk 



,2 



2 4M£ . r , r/ 
m - sm{M(t - t Q )} 



Defining the varying frequency as 



Xk^ . (22) 



k 2 



m 2 — 



4M£ 

(t - to) 



sm{M(t-t Q )} , (23) 



the same equation can be written as that of a parametric 
oscillator j(k + u\xk — 0. 

The particle production is achieved via the oscillating 
term in the above expression, through the already men- 
tioned parametric resonance. Inserting a variable z con- 
structed from the relation M(t — t Q ) = 2z ± ir/2, where 
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the different signs correspond to the sign of £, leads to 
the Mathieu equation 



d 2 Xfc 
dz 2 



+ [A k -2qcos(2z)} X k^0 , 



(24) 



where the parameters Ak and q determine the strength 
of parametric resonance and are given by 



4fc 2 



Am 2 



a 2 M 2 M 2 



M(t-t ) 



(25) 



The amount of parametric resonance can be described by 
a stability-instability map of the Mathieu equation [12l . 
[2 H3| , which have instability bands where the pertur- 
bations grow exponentially with different rates [2(| Hl| . 

In an expanding universe both Ak and q vary in 
time. Since the field passes many instability and stabil- 
ity bands, the growth changes with the cosmic expansion. 
The non-adiabaticity of the change of the frequency uik 



can be estimated by the quantity r n 
small k and to, becomes 



which, for 



B. Non-minimal coupling 

One addresses here a model that exhibits a non- 
minimal coupling between geometry and matter, as ex- 
pressed in the action Eq. (p}. Varying with respect to 
the metric leads to the modified field equations 



(27) 



{F 1 - 2F 2 p) ( iV - - 3ilv R ) = f 2 T^ 



-(/i - F^g^ + F 2 pRg^ + (Ft - 2F 2 p) , 

with Fi = fi(R); inserting the Robertson- Walker metric 
Eq. the temporal and spatial components read 



3(F 1 -2F 2 p)H 2 = (l + f 2 )p 

-\[h- (Fi - 2F 2 p) R] - 3Hd (Fx - 2F 2 p) 



(28) 



and 



M 



\cos{M(t-t )}/(t-t )\ 
|4A/£sin{M(t-t )}/(i-t )|3/2 



(26) 



The non-adiabatic regime corresponds to r na > 1. How- 
ever, as one can see, the condition M(t—t Q ) = mr leads to 
a more efficient non-adiabatic particle production, since 
fna 3* lj an d this occurs when the Ricci scalar vanishes. 

For instance, in the Starobinsky model @ the mass- 
less x particles are resonantly amplified for |£| > 3 [l2j], 
and the massive ones with to ~ M can be created for 
|£| > 10. These computations can be further extended 
to take into account non-linear effects, like the coupling 
between different modes, which can be important at the 
end of the preheating stage [U, [z| ; the backreaction of 
the Ricci scalar that emerges from the additional con- 
tribution of the energy density of the created particles 
to the background cosmological dynamics via the Fried- 
mann equation can also be considered, but with no sig- 
nificant impact on the reheating process [13 • 

One remarks that preheating is not a mandatory in- 
gredient of inflationary models: indeed, in the standard 
reheating scenario with £ = 0, particle production is still 
possible, albeit it is achieved via a perturbative decay of 
the inflaton degree of freedom (vis-a-vis the curvature, as 
discussed in the previous section and Ref. [23|]) into two 
X quanta 0, EH, H3 . 

From Eq. (|25|) . one sees that £ = implies q = 0, so 
that Eq. (|24|) describes a parametric oscillator driven by 
the expansion of the Universe alone; hence, the aforemen- 
tioned crossing of the broad parametric ressonance bands 
does not occur, rendering standard reheating rather inef- 
ficient when compared with preheating (see Ref. pB| for 
a recent discussion of the dynamics of thcrmalization). 



-2(F 1 -2F 2 p)H=(l + f 2 )(p + p)+ (29) 
(d d -Hd )(F 1 -2F 2 p) . 

As expected, GR is recovered by setting fi(R) = 2nR 
and f 2 (R) = 1. 

Resorting to the Bianchi identities, one concludes that 
the energy-momentum tensor may not be conserved in a 
covariant way 



h 



(30) 



In fact, as expected in the absence of the coupling, 
f'2(R) = 1) the covariant energy- momentum conservation 
is recovered. This feature implies that the motion of the 
matter distribution described by a Lagrangian density C 
does not follow a geodesic curve. Clearly, a violation of 
the Equivalence Principle may emerge if the r.h.s. of the 
last equation varies significantly for different matter dis- 
tributions, which suggests a method of testing the model 
and imposing constraints to the associated coupling con- 
stants. Of course, this putative violation is detectable 
only in astrophysical and/or cosmological contexts. 



III. INFLATIONARY REGIME 

A. Non-conservation of energy-momentum tensor 

In order to extend the standard preheating scenario in 
the context of Starobinsky inflation, one assumes a linear 
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coupling between matter and geometry 2 , as well as the 
usual Starobinsky curvature term, 



MR) 
MR) 



2k R 



R 2 
6M 5 



(31) 



1 + 2C 



R 

w 



where £ is a dimcnsionless parameter of the model and 
M has dimensions of mass. 

One considers that matter is described by a perfect 
fluid with the appropriate energy-momentum tensor, 



T, 



9 

p) U^U V + pg^ 




(32) 



where is the four- velocity, p is the density and p is the 
pressure, which is related with the former by the equation 
of state (EOS) p = up. 

Three different matter species are assumed, for com- 
pleteness: radiation, with L r = p r and u> r = 1/3, matter 
(taken as an ultrarelativistic particlc-antiparticle plasma) 
with C, n = —p m [26[ and ui m = 1/3 and a scalar field x 
with 



: Wx - V( X ) 



(33) 



In order to obtain the Friedmann equation, one re- 
quires the computation of d, obtained from the time 
component of the non-conservation law for the energy- 
momentum tensor of each species, as depicted below. 

Defining the energy density p x = —d^x^x/^ + V(x) 
and pressure p x = C x of the scalar field \i one gets a 
formally equivalent result for this component and for ra- 
diation, 



pi = - ( 3H + j; R ) (po+pj) = 

p r = -(4F+||)p r 
X = -(3H+%)x-V'(x) 



(34) 



(35) 



assuming that the curvature is high enough so that R > 
M 2 /(2£) -> f 2 (R) « 2£R/M 2 . 



2 Although not discussed here, it can be hinted from Eq. i51l that 
higher order additions f2(R) = 1 + £(R/M 2 ) n , with n > 1 would 
(in the fa(R) ~ 1 regime) deviate the dynamics from the desired 
parametric resonance — although a more convoluted preheating 
may perhaps still be achieved. 



As has been seen in the context of the accelerated ex- 
pansion of the Universe @ , one finds that the relativistic 
matter component is covariantly conserved, 



F2 

h 



T™ - g^C m ) V^R 



(36) 



-f (Too + 900Pm) R = 



and hence, 



Pm(t) 



-ZH(p n 

(CLi 
Pmi — 



(37) 



In the case of radiation, one may integrate directly the 
above expression to obtain the exact result 



p r (t) 





\MRiY 


4/3 


(a(t)) 


I MR). 





(a(f)) 


^M 2 ^ 


-2^1 




_M 2 - 


\-2iR 



4/3 



(38) 

4/3 



where p ri , a t and Ri are initial values of the radiation 
density, scale factor and curvature, respectively, and the 
final step assumes that f2(R) ^> 1 R 3> Af 2 /2£; con- 
versely, in the weak coupling regime /2OR) ^ 1 — > i? <C 
M 2 /2£ the above equation naturally reduces to the con- 
servation law p oc a~ 4 . 

Although both the matter field x an d radiation ad- 
mit the same Lagrangian density, C = p, the former is 
not characterized by a constant EOS parameter u> x = 
p(x)/p(x) — except for massless, non-self- interacting 
particles, V(x) = 0, where p x = p x — x 2 /2 —> uj x = 1, 
such that the scalar field behaves as ultra-stiff matter. 
One may compute the evolution of p x in the latter case, 
following Eq. (|3"5)) for x- 



Px( f ) = Pxi\T77\ 



a 



MRi) 



Pxi 



(t)' 



MR) 

Ri 
R 



(39) 



for an initial value of the scalar field density p X i. 

The dissimilitude between Eqs. (|3T|) and (f3"8")l is strik- 
ing, since in GR one naturally finds that matter follows 
the same conservation law as radiation, as they share the 
common EOS p = p/3. However, one cannot extend this 
common behaviour into the model here considered, i.e. 
assume that both components will obey the same non- 
conservation law. 

Indeed, one of the major features of the model embod- 
ied in action Eq. (HJ is that the Lagrangean density itself 
acquires a direct physical significance, which appears ex- 
plicitly in the modified field Eqs. (|2"T|) and the covariant 
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non-conservation law Eq. ([30)) — instead of only implic- 
itly, via the definition of the energy-momentum tensor. 
Hence, although the EOS of matter and radiation is the 
same, they follow distinct temporal evolutions because of 
the different Lagrangian densities, C r = p and C m = p. 

This said, one remarks that the choice of Lagrangian 
density for matter is the subject of some debate in the 
literature, mostly because it is never directly involved 
in computations in GR (see Ref. [26| and references 
therein). However, this discussion is not crucial to the 
results here presented, nor is the prescription C m = p 
paramount to their derivation: if one chose C m = C r = p 
instead, the non-conservation laws for matter and radia- 
tion would obviously be the same — but this would only 
alter the numerical coefficients affecting the r.h.s. of Eq. 
(|40p . not the ensuing discussion concerning the validity 
of the Starobinsky regime; from a less technical stand- 
point, one could invoke to dominance of radiation in the 
primitive Universe to ascertain the lack of importance of 
the choice for C m in the present context. 

Given the condition f 2 (R) > 1 -> R > M 2 /2£, 
one sees that the value of the "mass-scaling" parame- 
ter £ does not play a role in the above non-conservation 
laws. Furthermore, since the slow-roll regime yields 
R = —24H 4 e, with e <C 1, one does not expect the 
non-trivial contributions to make a significant impact, 
so that the three matter species evolve in a classical way. 
Finally, given the purpose of generalizing the preheating 
mechanism in the context of the non-minimal coupling 
model, one assumes that £ > 1, so that the condition 
R 3> M 2 /(2£) is always fulfilled during slow-roll (which 
ends when R ~ M ), and so the above non-conservation 
laws hold throughout this phase. 



where p = p r + p m + Px i s the total energy density. 

One assumes that the scalar field x has negligible ini- 
tial potential V(x) and kinetic energy x 2 , so that one 
may disregard it in the above equation; also, one assumes 
that the values for these quantities at the end of the 
slow-roll regime are still well below the dominant con- 
tribution, i.e. the scalar field will only play a relevant 
role in the (p)reheating phase. Recalling that R > M 2 
in the Starobinsky regime, the remaining conditions for 
Starobinsky inflation amount to 

^1 < M 2 <R , (41) 

K 

with j = r, m. Since ultrarelativistic matter behaves as 
radiation, one may use the Stefan-Boltzmann law pj oc 
gjT* (where gj is the number of degrees of freedom and 
T is the temperature, assuming thermal equilibrium) to 
write 

£«<5gWa;" ,ir ' rl/ '- < 42) 

after using that Mw3x 10~ 6 Mp @. 

As stated before, one assumes that £ > 1, so that this 
parameter does not play a role in the evolution of matter 
during slow-roll, but it becomes relevant in preheating, 
at the onset of the condition R ~ M 2 /(2£). This, cou- 
pled with the expected smallness of the R/R terms in Eq. 
([33]) , leads one to consider that the temporal evolution of 
the density of the considered matter species does not dif- 
fer significantly from that expected in standard Starobin- 
sky inflation. Indeed, using the slow-roll approximation 



B. Starobinsky regime 



R = 6H 2 (2 - e) -> R ~ -24£T 3 e 



(43) 



One aims to obtain a regime of Starobinsky inflation, 
where the terms derived from the curvature contribu- 
tion /i(i?) = 2k(R + R 2 /{6M 2 )) dominate, thus lead- 
ing to the inflationary evolution described in the previ- 
ous section. For this, one should rederive the Fricdmann 
equation, taking into consideration also the effect of the 
non-minimal coupling f2(R) = 1 + 2£R/M 2 , and set the 
conditions which should be fulfilled so that these are sub- 
dominant with respect to the latter. 

This said, using the — component of the Ein- 
stein equations and the (non-) conservation laws obtained 
above, one gets the modified Friedmann equation, 



2g 

iM 2 



MI 2 



3p„ 

1 

M 2 

kR 
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R 2£ /4 



■X-v 



(40) 



6 + T 3^ 



2 X 



H 



3 Pr 



- 1 



HR+f- 

OK 



one gets 



so that 



R 



-2HRe 



Pr 



(3 - 2e) Hp 



X = -(3-2e)Hx-V'( X ) 



(44) 



(45) 
(46) 



and, since e < 1, the contribution from the non-minimal 
coupling may be disregarded, and the energy-momentum 
tensor of all matter species is almost completely con- 
served (see paragraph MI Cl for more detail). In partic- 
ular, this means that no decreasing effect on the cool- 
ing rate is attained, i.e. the Hubble parameter always 
dominates the r.h.s. of the considered equation and the 
density (and temperature) never rises due to the non- 
minimal coupling. Hence, it suffices to evaluate Eqs. (|41[) 
and (|42[) at their maximum initial values. 
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Upper bounds for £ can be obtained by assuming that 
the initial inflationary temperature is, at least, of the 
grand unified theory (GUT) scale T GUT ~ 10 15 GeV ~ 
10~ 4 Mp. Hence, one obtains 

{ < 10-« (J^V ~ 10* . (47) 

The obtained range 1 < £ < 10 4 can be translated into 
an equivalent mass scale appearing in the non-minimal 
coupling term f 2 (R) = 1 + R/M 2 , with M 2 = ~ 
10~ 2 A/-M; hence, one does not introduce a strong mass 
hierarchy in the non-trivial components of the models, 
and may consider only the mass scale M as the key in- 
gredient to the model. 

Finally, one sees from Eq. (|42|l that the maximum 
initial temperature compatible with £ > 1 is given by 
T m ax ~ 10~ 3 Mp ~ 10 Tqut', hence, our model seems 
to be particularly well suited for a GUT scale inflation, 
although it fails for Planck scale inflation. The latter, 
in turn, is not quite desirable on account of most of the 
baryogenesis mechanisms (see however, Ref. p7j|). 

The above constraint is not regarded as a caveat of the 
present approach, since it is not in contradiction with 
the required minimum number of e- folds N > 60: in- 
deed, inserting Hi ~ T ~ Tqut into the result from 
the Starobinsky regime N ~ 3(H/M) 2 , one obtains 
N = 10 3 - 10 5 > 60, having again used the CMB value 
M ps 3 x 10~ 6 M P . 



C. Slow-cooling due to the non-minimal coupling 

As argued before, it is not feasible to use the effect of 
the non-minimal coupling to alleviate the supercooling 
ensued after inflation. The impossibility of attaining this 
"slow-cooling" regime is assessed by simply writing the 
final radiation density after slow-roll, as obtained from 
the direct analytical solution Eq. (|38[l : 

^/ = ^(^) 4 (^) 4/3 = P ri e-^(4^/3 i (48) 

having used R f = 6Hj = M 2 and R t ~ 12H? ~ 
ANM 2 = 4NRf. Hence, one may only obtain an in- 
crease in density by a factor (4iV) 4 / 3 (when compared 
with the value obtained without the effect of the non- 
minimal coupling) or, using the Stefan-Boltzmann law 
p r oc T 4 , a factor (4./V) 1 / 3 in temperature. Consider- 
ing that the number of e-folds is constrained to be of 
the order 10 3 — 10 5 , due to the validity conditions of 
the Starobinsky regime already discussed, one obtains a 
compensating factor of order 10 — 100, manifestly insuffi- 
cient to end inflation with a "warm" Universe. Thus, one 
must deal with the issue of implementing reheating after 
slow-roll, as will be considered in the following section. 



IV. REHEATING 
A. Starobinsky evolution 

Assuming that the condition Eq. (|4Tj) holds, one may 
approximate the modified Friedmann equation (|40[) by 
the one arising in standard Starobinsky inflation, so that 
the slow-roll treatment gives rise to the temporal evolu- 
tion already discussed. 

Given that £ > 1 , the end of slow-roll e ~ 1 corresponds 
to R ~ M 2 > M 2 /(2£); hence, one might assume that 
the non-minimal coupling f 2 (R) = 1 + 2£i?/M 2 would 
take hold of the dynamics until it too gets close to unity, 
when R < M 2 /(2£). However, one should keep in mind 
that its effect is mediated by the value of the density 
of the considered matter species, according to Eq. (|41[) : 
since this condition stems from the modified Friedmann 
equation, Eq. (14TJ1) . regardless of the validity of the slow- 
roll condition e <C 1, and that the density has dropped 
exponentially during inflation. Hence, one concludes that 
the post-inflationary dynamics is achieved by disregard- 
ing the terms due to the non-minimal coupling: the evo- 
lution of the Hubble parameter and the curvature is then 
given by Starobinsky model post-inflationary phase de- 
scribed by Eqs. dTTJ) and (T2"|l. 

B. f-2 ~ 1 regime 

As stated before, one attempts now to use the dynam- 
ics of the matter scalar field x to provide the required 
reheating of the Universe, which after inflation is essen- 
tially at zero temperature. For concreteness, one dubs 
the timespan during which f 2 (R) ps 2£R/M 2 as the cou- 
pled regime, while f 2 (R) ps 1 is naturally termed as un- 
coupled. 

Clearly, since the curvature oscillates, one cannot as- 
sume that the coupled regime prevails at the onset of os- 
cillation, as during each of these the curvature vanishes 
and one has f2(R) = 1- Furthermore, since its amplitude 
also drops with t _1 , at t = t Q + 2£/M (that is, after £/-7T 
oscillations) it is below M 2 = M 2 /2£, a point at which 
one can no longer assume that f 2 {R) ps R/M$ ^> 1. 

From Eq. (|4"T|) , the number of oscillations before one 
crosses over to the f 2 (R) ps 1 regime is constrained to 
be less than ~ 10 3 ; if £ is close to unity (and no hierar- 
chy exists between M and M 2 ), not even one oscillation 
is completed. Bearing this in mind, one considers only 
this uncoupled regime; this can be considered a worst 
case scenario, given that one neglects any reheating that 
might occur during the coupled regime. A numerical 
analysis (not performed here) should allow for a better 
discrimination of both regimes and assess the impact on 
the dynamics and reheating. 

However, disregarding this "coupled reheating" should 
bear no relevance on the final result, since the reheating 
generated during the uncoupled regime should be much 
more important: indeed, this arises from the term F 2 / f 2 
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appearing in Eq. (|35|) . which (since F2 = const.) rises as 
the curvature falls and f2(R) drops to unity. 

As a result of the above discussion, one proceeds to 
consider the fz{R) ~ 1 uncoupled late-time regime as 
the only channel to implement a suitable reheating mech- 
anism in the model under scrutiny. Naturally, this does 
not coincide with the minimally coupled regime f2(R) = 
1, since the kinetic term F2 is non- vanishing — only the 
covariant non-conservation laws for matter and radiation 
collapse to their GR counterparts, as Eq. ([37J) becomes 
identical to Eq. (|38| . i.e. p m oc p r oc a~ 4 . 

One begins by noticing that Eq. (f35|) no longer applies, 
since -F2//2 —> 2£/M 2 . As before, one may resort to the 
analytical integration of the non-conservation law for the 
radiation pressure or a massless scalar field with constant 
EOS parameter w x = 1, Eq. ((38)) : 



pj(t) = p 2 j 
Poj 



a(t)J 

o(t 2 ) \ 3(1+ " j) 
a(t) 



HR2) 



MR) 



(49) 



where t — t^ refers to a moment at the onset of the regime 
R < M 2 /(2£) -> f 2 (R) « 1. One obtains the standard 
GR density evolution of a fluid with constant EOS pa- 
rameter u>j, and thus one concludes that no reheating 
occurs due to the non-minimal coupling between curva- 
ture and matter (which is always covariantly conserved) , 
irrespectively whether radiation or a massless scalar field. 



C. Preheating 

One is left with the possibility of a non- vanishing po- 
tential V(x)', one assumes that V(x) = m 2 x 2 , so that 
X is a quantum field which may be decomposed (in the 
Heisenberg representation) as in Eq. (fl"9]l . such that each 
Fourier mode Xk obeys the differential equation 



Xk 



3H 



Fo 



R)Xk 



„2 



Xk =0 



(50) 



Introducing the field redefinition Xk = a 3 / 2 f^^Xk an d 
the new variable 2z = M(t — t )±ir/2 (depending on the 
sign of £), one can absorb the x'k term, obtaining 



IF2 
2/2 



/ 2k 
[aM 
IF2 
2 ft' 



2m\ 2 


H' 


H 2 


mJ 




HR' 
6 M 


-R"^j 


x k = o 



(51) 



with the prime denoting differentiation with respect to 
z. Notice that the redefined field X k closely approaches 
the one obtained via the transformation performed in 



standard preheating scenarios, Xk Xok = a 3 ^ 2 Xfe as 
the curvature drops and /^(-R) — > 1- 

Moreover, one could choose to perform the latter trans- 
formation, thus absorbing only the term 3Hxk- The re- 
sult would have yielded a similar equation to Eq. (|51[) . 
but with an added x'k term involving the non-minimal 
coupling factor i^/Za- A similar consideration was put 
forward in another study, in the context of reheating via 
a matter scalar field with a non-standard kinetic term 
0. 

The generalization is obvious: while standard preheat- 
ing resorts to a non-minimal coupling in the mass term 
only (which enables an oscillating effective mass term 
We// — m2 + £R) an d the latter study considers a non- 
minimal coupling to the kinetic term only, our proposal 
does not distinguish between the components of the La- 
grangcan density (or, in fact, between Cj for different 
matter species). 

In the fi{R) ~ 1 regime, Eq. (|5l"j) reads 



x k 



( 2k Y 

M 2 VM 2 



2rri 
HR' 



H' _ H 2 



(52) 



6- 



M 



R" 



Xi 



In the late-time approximation z» 1, one can write 



ff(z)~^[l + sin(2z)] 
iz 



(53) 



R{z) ~ 3MH'(z) 



2M< 



cos(2z) , 



so that the terms due to the Hubble expansion inside the 
brackets of Eq. (f5Tj) read 



3 — = -cos(2z) , 

M z v ' 



(54) 



M 2 



1 + sin(2z) 



and H' ^> H 2 /M, as expected. The terms arising from 
the non-minimal coupling become 



M 2 ) 
HR' 



M 3 
R" 



l^sin 2 (2z) , 

z z 

-^|sin(2z)[l + sin(2z)] , 



(55) 



z 
8i 



M 2 



cos(2z) 



Since the uncoupled regime implies R < A/ 2 /(2£) — > 
z/(2£) > 1, one concludes that the last term £R"/M 2 
dominates the other two arising from the non-minimal 
coupling (although different in scope, a discussion con- 
cerning the composite effect of different harmonics in the 
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Mathieu Eq. ([53]) can be found in Rcf. [LJ). Further- 
more, one sees that the latter compares with the domi- 
nant term from the Hubble expansion as 



W 3_H> 

q M 2 5 M 



(56) 



indicating that the non-minimal coupling dominates the 
dynamics for £ > 1. Hence, Eq. (|52]l simplifies to 



M 



•M 2 



X k = , (57) 



which, inserting Eq. (|55|). can be cast as a Mathieu equa- 
tion, 



[A k -2qcos{2z)]X k = 



(58) 



with 



A k = 



f 2k\ (2m 



4£ 

q = — 

z 



(59) 



This is exactly the same equation as in the usual pre- 
heating scenario obtained via an ad-hoc coupling term of 
the form £i?x 2 ■ Thus, one may extrapolate directly from 
the results already available for this reheating model. In 
particular, one recalls that massless particles may be pro- 
duced for values of the coupling parameter as low as 
£ > 3, and massive particles for £ > 10. Hence, the 
allowed range 1 < £ < 10 4 enables an ample production 
of massive particles. 



V. DISCUSSION AND CONCLUSIONS 

In this work one demonstrates that a generalized cou- 
pling between matter and geometry naturally extends the 
preheating mechanism in the context of the Starobin- 
sky inflationary model. This novel preheating scenar- 
ios arises from a scalar field x with a scalar curvature- 



dependent effective mass term tu 2 ^ 



£i? or 



with a non-canonical kinetic term. One finds that a 
generalized coupling of the form f2(R)£ (with a lin- 
ear f2(R) = 1 + 2£i?/M 2 ) leads to a dominant term of 
the form R/M 2 in the Mathieu equation that describes 
the quantum production of particles with momentum fc; 
given that R(t) is rapidly oscillating in the reheating era, 
one gets R/M 2 ~ R, and the aforementioned preheating 
mechanism is fully recovered. 

The generalized coupling implies that one could also 
consider a possible dominance of the terms arising from 
/2 in the modified Friedmann equation, since these are 
coupled not only with the pressure and density of the 
scalar field x^ but also with matter and radiation. One 



finds that constraining the dynamics so that a Starobin- 
sky inflationary era occurs translates into a mild condi- 
tion for the coupling parameter 1 < £ < 10 4 , more than 
sufficient to allow for an efficient production of massive 
particles, thus leading to a successful reheating of the 
Universe. 

The obtained result should not be viewed as a sole ex- 
ample of the explanatory power of the generalized non- 
minimal coupling model under scrutiny, but regarded to- 
gether with other available results: namely, the descrip- 
tion of the galaxy rotation curves via a mimicking "dark 
matter" mechanism Q, and the possibility of account- 
ing for the current accelerated expansion of the Universe 
Q. As discussed in those works, the lower values of R 
at the late-time and large-distance regimes where these 
phenomena arise, respectively, naturally lead to the con- 
sideration of an inverse power-law coupling f2{R) ~ R~ n , 
with n > 0. 

On the contrary, considerations regarding inflation 
must be implemented via a direct power-law coupling (in 
this work, a linear one), which guarantee that its effect 
does not grow after reheating, leading the Universe away 
from the desired radiation-dominated era. 

Conversely, the inverse power-laws considered in Refs. 
[1, Q bear no significance in the present context, given 
that the curvature was sufficiently high during the infla- 
tionary era to render its effects completely negligible. 

In order to unify these different couplings, one should 
consider that each regime (early vs. late time, central vs. 
long range) is dominated by the influence of a particu- 
lar contribution to an overall, yet unknown non-minimal 
coupling function, which can then be written as a Lau- 
rent scries, 



f2(R) 



2 (m= 



(60) 



Naturally, other particular phenomena and environ- 
ments may be of interest, if the typical values of the 
curvature are such that yet unprobed terms of the above 
series dominate the dynamics. 

Finally, one discusses another context where the influ- 
ence of a linear coupling f2(R) may be addressed: the 
central regions of a spherical body (such as a star), as 
its core might be sufficiently dense to constraint its effect 

M- 

Using the hydrostatic equilibrium of the Sun as a test 
case, it was found that the central density and available 
observational precision of central observablcs yield only 
a modest bound for the non-minimal coupling strength; 
translated into the notation of this work, one gets £ -c 
10 78 , several orders of magnitude above the upper limit 
of the range obtained here, 1 < £ < 10 4 . 

With the above discussion in mind, the current work 
may be regarded as yet another effort to reconstruct the 
full form (|60l) of the non-minimal coupling — in this case 
by constraining its linear term in a far more effective 
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fashion than through the aforementioned study of solar 
hydrostatic equilibrium. 



Appendix A: Multi-scalar field formulation 

It is widely known that f(R) theories can be rewrit- 
ten as the Einsteinean gravity plus a scalar field model, 
with the curvature being dynamicall y id entified with a 
scalar <f> = R (or a function thereof) [23|. Similarly, the 
discussed non-minimally coupled model Eq. ([1]) is analo- 
gous to a multi-scalar field theory, with two scalar fields, 
<f> = R and if) = £, as is illustrated below. One may per- 
form a conformal transformation to the Einstein frame, 
so that the curvature term appears decoupled from these 
fields, and redefine the latter, obtaining the equivalent 
action 



S = J V^gd 4 x x (Al) 
2k [R ~ 2g^a ij ^ v - iU(<p\ p 2 )] + f 2 (p 2 )C* 
where p 1 and p 2 are scalar fields, er^ is the field-metric 



1 1 

-1 



the potential is given by 




(A2) 



(A3) 



AO 2 

2k 



■ cxp 



2\/3 ! 



and C* = cxp[-(4 % /3/3)^ 1 ]. 

The two scalar fields are related with the scalar curva- 
ture and the non-trivial fi{R) and f2(R) functions as 



i \/3. 
f = ~2 lo S 



F 1 (R)+F 2 (R)C 



2k 



ip 2 = R . (A4) 



Inserting the expressions fi(R) = 2k (i? + R 2 /(6M 2 )) 
and f 2 (R) = 1 + 2£R/M 2 , one gets 



Vs. 

T log 



3M 2 



(A5) 



and the potential 



Appendix B: Evolution with trivial fi(R) term 

In this section, one attempts to obtain an inflationary 
solution to the modified Friedmann Eq. (|40|) by resort- 
ing only to the non- minimal coupling f 2 (R), setting 
fi(R) = 2kR, its GR form, and taking a general non- 
minimal coupling f 2 (R). Considering only the dominant 
matter contribution of radiation, p w p r , for simplicity, 
the Friedmann equation then becomes 



l-Fo 



' 3k 
F 2 R 



H 



(Bl) 



P_ 
6/t 



4f; 



J1 + ZzJL - fL I -^-HR 



3 /• 



3k 



If one assumes that the effect of the non-minimal cou- 
pling dominates the usual GR terms, the above can be 
approximated by 



1 / F 2 R 
2\ h + — 



HR . (B2) 



Remarkably, in this strongly coupled regime the den- 
sity does not appear in the dynamical description of the 
cosmological evolution, and the obtained Friedmann Eq. 
is strikingly similar to the analogous expression arising 
from usual f(R) theories This observation paves the 
way to the possible use of the non-minimal coupling f 2 
as the driving force behind inflation — instead of serving 
only to implement a reheating mechanism, with the infla- 
tionary regime occurring due to the fi(R) contribution. 



1. Power-law non-minimal coupling 

With the above discussion in mind, one first considers 
couplings of the form 



( ) 

F 2 (R) « ^h(R) , 



(B3) 



assuming that the characteristic mass scale M 2 is well 
below the initial value for the scalar curvature, R 3> M 2 - 
Clearly, the exponent n must be positive, so that a de- 
creasing curvature will eventually make the non-minimal 
coupling terms subdominant with respect to the GR con- 
tributions. 

Substituting into Eq. (|B2|) . one gets 



[/(^ 2 ) = ^ 2 exp(-^V) x (A6) 



6M ; 



cxp 



2%/3 l 



H 




(B4) 



Notice that the mass scale M 2 vanishes from the above; 
resorting to the slow- roll approximation R = 6H 2 (2 — 
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e) —> R ~ — 24H 3 e, one may promptly obtain the slow- 
roll parameter 



36 + 18n 
36 + Tin + An 2 



(B5) 



The linear coupling n = 1 yields ei = 0.80, which (al- 
though smaller than unity) does not fulfill the slow-roll 
condition e n <C 1. Indeed, only for the unnatural range 
?i > 40 does one get e n < 0.1. Also, notice that in the 
large n regime one obtains = 9/2n. 

One now checks that the evolution of the radiation den- 
sity p; using Eq. (|35[) and the slow-roll approximation, 
one has 



(B6) 




Hp = 



having used Eq. (|B5|) . Only for n < 2.32 does one gets 
a diluting Universe, p < 0; it is trivial to check that, for 
large n, p = 8Hp — > p oc a 8 . Hence, the slow- roll range 
?i > 40 yields an exploding density, which would push 
back the expanding Universe into Planck scale concen- 
tration, rendering the model inapplicable. 

Even if the density does not become radically high, 
another problem would arise. Firstly, it can be shown 
that, even though the density increases, the f 2 (R) de- 
rived terms in the modified Fricdmann Eq. (|Blj) di- 
minish: this occurs since the latter are of the generic 
form F 2 Rp, and both F 2 and R decrease. If this was 
not the case, then the Universe would become eternally 
dominated by the non-minimal coupling, even after the 
slow-roll condition e<l fails — possibly including the 
so-called super-inflation, H > 0. 

However, this would only allow for inflation to end and 
the standard GR evolution to ensue: but this exit would 
be anything but graceful, as the Universe would be ex- 
tremely hot, instead of at essentially zero temperature, 
overshooting the reheating temperature by several orders 
of magnitude. Hence, one concludes that no power-law 
models can solely drive a physically meaningful inflation- 
ary regime. 



2. Exponential non-minimal coupling 



Following the discussion after Eq. (|B2[) . one now re- 
sorts to an exponential non-minimal coupling function, 



MR) 



exp 



R 

Ml 



> 1 



(B7) 



again considering the coupled regime R 3> Mf . This is a 
somewhat natural choice, since it collapses to the trivial 
GR scenario for small values of C, /2(A) — > 1; as before, 
one considers a positive exponent n. Substituting into 
Eq. (1B"2]) . one gets 



2 



R 



(B8) 



R 

Ml 



1-n 



HR 



Considering that R ^> Mf , the above simplifies to 



R 
6" 



R 

Ml 



HR 



(B9) 



One uses the expansion 



R 

Ml 



6H 2 (2 



Ml 



(B10) 



V M. 



which, together with the usual slow-roll approximation, 
yields the slow-roll parameter 



_9_ / M 2 \ 
2n \2V3Hj 



2ii 



(Bll) 



having again considered that H 3> M 2 . Given a suffi- 
ciently high initial Hubble parameter Hi, the the slow-roll 
condition e„ <gc 1 is clearly valid for a positive exponent 



Inspection of Eq. (|B11|) shows that the slow-roll phase 
ends when 



l/2n 



H = f n M 2 , /„ = 2V3 I- , (B12) 

with the numerical factor satisfying /„ = 2\/3, for n — > 
00. 

A striking result is that, if one takes n = 1, i.e. a 
simple exponential coupling, then 



3 (M 2 



(B13) 



which, from Eq. ([8} , is identical to Starobinsky inflation, 
redefining M 2 = 2M/3. 

As in the previous paragraph, one now ascertains what 
is the evolution of the radiation density p, following Eq. 



12 



(|35|) . As it turns out, the obtained result does not depend 
on the exponent n, as seen below: 



4F 2 



(B14) 



H 



R 



R 
R 



■in 



2n ( 2V3H \ ' 



Hp = 8Hp > 



Again, one encounters the same pathology already dis- 
cussed: the density increases to pre-inflationary values; 
even if inflation ends before it reenters the Planck scale, 
one would be left with a burning hot Universe, much 
higher than the required reheating temperature. 

As a final remark, one concludes that a suitable infla- 
tionary regime leading to a final reheating temperature 
appears unattainable resorting solely to a non-minimal 



coupling, for a wide choice of "natural" forms f2(R) (not 
shown in this manuscript), unless more model parame- 
ters are considered — which would have to be introduced 
in an ad-hoc and unphysical fashion, in order to alleviate 
the above issue and tame the temperature evolution of 
the Universe. 

Nevertheless, the mimicking behaviour exhibited in 
Eq. (|B1|) . attested by the recovery of Starobinsky infla- 
tion in the /^(-R) = exp(iJ/Mf) case, could perhaps be 
useful in the context of the early Universe, although none 
is envisaged in this study; a similar mimicking behaviour 
was obtained when considering the current accelerated 
expansion phase Q . 

The attempts described in this Appendix, although un- 
successful, serve to better frame the main purpose of this 
study: to show that a generalized non-minimal coupling 
may be paired to the Starobinsky quadratic curvature 
term in a modified Einstcin-Hilbcrt action, so that the 
former successfully enables the preheating of the Uni- 
verse, while it abstains from influencing the inflationary 
dynamics — which is driven solely by the latter. 
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